The dynamics of apparent horizons in Robinson- Trautman spacetimes 
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We present an alternative scheme ol finding apparent horizons based on spectral methods applied 
to Robinson- Trautman spacetimes. We have considered distinct initial data such as representing 
the spheroids of matter and the head-on collision of two non-rotating black holes. The evolution of 
the apparent horizon is presented. We have obtained in some cases a mass gap between the final 
Bondi and apparent horizon masses, whose implications were briefly commented in the light of the 
thermodynamics of black holes. 
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I. INTRODUCTION 



One of the most important problems in classical 
General Relativity is the evolutijjn of apparent hori- 
zons. The apparent horizon [l| is defined as the 
outermost marginally trapped surface that can be 
located on each spacelike surface during the over- 
all dynamics of the spacetime. According to cosmic 
censorship hypothesis there must exist outside the 
apparent horizon an event horizon T], and for this 
reason apparent horizons are the key structures that 
signalize the formation of black holes in gravitational 
collapse, as well play relevant role in the merging of 
black holes Besides the apparent horizon an- 
other typical structure present in a spacetime that 
contains a black hole is the event horizon, but its de- 
termination depends on whole history of the space- 
time due to the fact that an event horizon is the 
boundary that separates those null geodesies that 
reach infinity from those that not. In essence, while 
the event horizon is a global structure, the apparent 
horizon is local meaning that it can be determined 
at each instant. Therefore, the construction of ap- 
parent horizon finders is a crucial issue in numerical 
relativity that has received a great deal of attention 
in the last years 0]. Basically, these codes are built 
to solve numerically the apparent horizon equation, 
which is a nonlinear elliptical equation, simultane- 
ously with the numerical evolution of the spacetime. 

In general most of the numerical strategies to solve 
the apparent horizon equation are based on the fi- 
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nite difference techniques. On the other hand, nu- 
merical codes based on spectral methods jl, Q have 
increased considerably in the last years mainly due 
to the economy of the computational resources to 
achieve a desired accuracy. In this direction we shall 
present here a simple and efficient numerical strategy 
using a convenient combination of Galerkin [6|] and 
collocation methods 0, [1] to determine the evolu- 
tion of the apparent horizon of Robinson- Trautman 
spacetimes 

The Robinson- Trautman (RT) spacetimes are the 
simplest class of asymptotically flat geometries ad- 
mitting gravitational waves with two interesting ba- 
sic features: (a) a RT spacetime can be interpreted 
as describing the exterior geometry of a bounded or 
isolated system emitting gravitational waves; (b) for 
regular initial data the RT spacetimes evolve asymp- 
totically towards the Schwarzschild black hole [l(j |. 
For the sake of completeness, the line element of the 
Robinson- Tratuman spacetimes can be conveniently 
written as 



ds^ = \{u,e) - 



2mo 



K \ 

2r— \ du^ + 2dudr - 
K / 



(1) 



where dot means derivative with respect to the null 
coordinate u, r is the radial coordinate, {9, ip) are the 
usual angular coordinates, and mg is an arbitrary 
constant. The Einstein equations can be cast in the 
following form 
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Here the subscript 9 denote derivative with respect 
to the angle 9; the function A(u, 6) is the Gaussian 
curvature of the surfaces (u = const., r = const.). 
The structure of the field equations is typical of a 
characteristic problem T7|, in which the first equa- 
tion is a hypersurface equation relating the functions 
X{u,9) and K{u,9), whereas the second equation is 
the evolution equation. Accordingly, once the initial 
data K{uo,9) is prescribed on a given null surface 
u — uq, the hypersurface equation fixes A(uo, 9), and 
the evolution equation determines K{u, 9) on the 
next null surface, and whole process repeats provid- 
ing the evolution of the spacetime. 

The only known analytical solutions of the 
RT field equations are the two forms of the 
Schwarzschild solution described by 



and therefore indicating a simple example of non- 
spherical collapse. However, these geometries do not 
have a future apparent horizon characterized by the 
vanishing of the null expansion associated to out- 
going future directed rays, but only past apparent 
horizon [3, [l^, [l^ . A past apparent horizon is the 
outermost boundary of past-trapped surfaces corre- 
sponding to that value of u; more precisely, consider 
a hypersurface S defined by 5 = r — V{u, 9) = 0, 
and in particular if 5 is a marginally past-trapped 
2-surface, the ingoing normal null vector — daS 
has vanishing divergence 



= 0. 



(7) 



K = Kq = constant. 



(4) 



From this equation it can be shown [1^ that the 
function V{u, 9) satisfies the following equation at 
each slice u — constant, 



K{9) 



cosh 7 + cos 9 sinh 7 



(5) 



The first reproduces the Schwarzschild black hole 
with mass Mbh = moKQ, and the second expres- 
sion a boosted black hole with constant velocity 
V = tanh7 with respect to an inertial observer at 
infinity. In this case the total mass-energy content 
is given by 



Mbh = 



Too if 



(6) 



Notice that this above expression corresponds to the 

total relativistic energy of a moving particle with 
— 3 

rest mass nioKo and velocity v. 

This paper is divided as follows. In Section 2 we 
present the apparent horizon equation and the nu- 
merical strategy based on spectral methods adopted 
to solve it. In Section 3 we exhibit the numerical re- 
sults that consists in testing the code along with the 
dynamics of the apparent horizon corresponding to 
initial data representing spheroids and the col- 
lision of black holes [ij] in RT spacetimes. Finally, 
the final remarks are presented in Section 4. 



II. SOLVING THE APPARENT HORIZON 
EQUATION USING SPECTRAL METHODS 

As already mentioned, RT spacetimes have inter- 
esting features such as the asymptotic flatness and 
the presence of gravitational waves, which can be in- 
terpreted as arising from a bounded distribution of 
matter evolving towards a Schwarzschild black hole, 



sin 9 



V 



2mo 



= 0, 



(8) 



which is known as the apparent horizon equation. 
The dynamics of the apparent horizon is obtained 
after solving this equation at each hypersurface u — 
constant, where the function K{u,9) is determined 
from the evolution equation There are few ana- 
lytical results about the properties of past apparent 
horizons in RT spacetimes. Tod has shown the 
validity of the isoperimetric inequality and the ex- 
istence of a unique marginally past-trapped surface 
at each hypersurface u=constant. 

The apparent horizon equation ([5]) will be solved 
here using a numerical scheme based on a suitable 
combination of Galerkin and collocation methods in 
a similar way we have implemented to solve the field 
equations ^ and ([3]). For this reason we shall briefly 
outline our previous numerical scheme p^l [isj for 
solving the field equations and, in the sequence, the 
procedure employed to integrate the apparent hori- 
zon equation. 

According to Ref. [3 the first step is to estabhsh 
the Galerkin expansion for the function K{u,9), 



Kliu,x) = e-^-t"'-) = exp ] , (9) 

\fc=o / 

where the subscript a indicates an approximation 
of the exact K{u,x). The angular coordinate 9 is 
replaced by a; = cos 9, N is the truncation order that 
indicates where the series stops, and the N+1 modes 
bk{u) are unknown functions of u to be determined; 
the Legendre polynomials Pk{x) were chosen as the 
basis or the trial functions. Next, an approximate 
expression for the function A(w, x) is obtained after 
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substituting the above expansion into the constraint 
equation or 

Xaiu, x) = e-«-("'^' + M^±il6,(„)P,(a;) j (10) 

These last two equations are substituted into Eq. 
to yield what is know as the residual equation 
associated to the evolution equation, 



and after performing the TV + 1 integrals, a set 
of ordinary differential equations for the modes hk [u) 
arises. Therefore, evolving these equations means to 
determine the dynamics of RT spacetimes since the 
function K{u, x) can be reconstructed at each u. 

The past horizon equation ([5]) will be solved at 
each time level u by applying a similar combina- 
tion of spectral methods as described previously. We 
have followed Ref. [15] and introduced an auxiliary 
function x) by setting 



N 



ResK{u,x) = 6mo ^ bk{u)Pk{x) 



V = 2mo exp(-F), 



(14) 



fc=0 

-Qa {U,X) 



in order to eliminate mo from the apparent horizon 
(11) equation. A natural Galerkin expansion for i^(M, x) 
is given by 



where prime denotes derivative with respect to x. 
Notice that the residual equation does not vanish 
exactly due to the adopted approximations for the 
functions K{u, x) and A(u, x), but as we have shown 
it converges to zero as the truncation order N is 
increased [l^. Following the Galerkin method, the 
projections of the residual equation with respect to 
a suitable set of test functions ^'n(x) vanish, namely 



M 



Fa{u, a;) = ^ ak{u)Pk{x), 



(15) 



k=0 



(Resi<-(u, a;), ^'„(x)) = / Resi<-(u, a;)vl>„(a:) da; = 0, 



(12) 



where M is the truncation order not necessarily the 
same as N (see Eq. (|9])). The apparent horizon 
equation is rewritten in function of F{u,x), and af- 
ter substituting the above expansion together with 
the approximate expressions for K{u, x) and A(u, a;) 
(Eqs. ^ and (fTU]) ). we have obtained the residual 
equation associated to the apparent horizon equa- 
tion 



for n = 0, 1, ..N. It means that the modes bj{u) are 
chosen in such a way that the residual equation is 
forced to be zero in an average sense [l^ . Following 
the Galerkin method we have selected, for the above 
integrations, the test functions to be same as the 
trial functions, ^'^(a;) — Pn{x). At this point we 
have introduced an additional approximation for the 
exponential term given by 



ResAH(^^,a;) = [(1 - a:^)^^^]' + 1 + ^ -fc(fc + 1) 



A;=0 



bu{u)Pk{x) - exp(Fa -I- Qa). 



(16) 



exp(-Qa(",a;)) 



N 

E 

.7=0 



(13) 



where Tj ix) is the Chebyshev polynomial of order j 
and N indicates the number of modes Cj . Basically, 
the motivation behind this approximation is to allow 
rapid and direct integrations of the residual equa- 
tion. As a consequence of the above expansion, the 
N + 1 modes Cj are related to the TV -|- 1 modes bk by 
assuming that the projections of the residual equa- 
tion ResQ(u,a;) = exp{-Qi,{u, x)) - J2f=o ^j^^ji^) 
with respect to the test functions 5'„(a;) = S(x — Xn) 
vanish, where xo,xi, .., a;^ are the collocation points 
associated to the Chebyshev polynomials. The ad- 
ditional approximation (jl3p is introduced into Ec^. 



As we have described before, the next step is to 
impose that all projections of the residual equation 
with respect to each basis function, Pn (x) , n = 
0, 1, M, must vanish. Schematically, we have 



(ResAH,-Pn(a;)) = J ResAHiu,x)Pnix) = 0. (17) 

Notice the presence of two exponential terms in the 
residual equation p6)) that can be reexpressed using 
additional approximations as. 



M 



exp(i^a('«, a;)) « ^ akTk{x), (18) 

k=Q 
M 

exp(Qa(w, x)) « PMx), (19) 



k=0 
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where ak and Pk are the modes associated to 
these new approximations, and M is the trun- 
cation order for both expansions. The pro- 
jections of the corresponding residual equations, 

ResF = exp{Fi,{u, x)) - J2kLo otkTk{x) and Resg = 

exp((3a(w, a;)) - Xifclo PkTkix), with respect to the 
test functions S{x — Xn), with Xn = 0, 1, .., M being 
the collocation points of Chebyshev polynomials, are 
forced to vanish. Consequently, two sets of M -I- 1 
algebraic equations relating the modes (a^, [3k) with 
(aj,bk), respectively, are generated. These approx- 
imate expressions are then inserted into Eq. (jl7p . 
yielding 



(ReSAH, Pr 



+ 1+-X 



AT M 

^ k{k + i)hk{u)Pk{x) - Mu)Pj{u)n{x) X 

Tj{x)}P„{x)^0. (20) 



After performing the above integrals, a set of M -I- 1 
algebraic equations of the type fk{aj,bj,ai, Pi) — 
is obtained. Since we can express the modes ak 
and Pk in terms of aj and 6j, and the modes bk are 
known at each m, we can solve, in principle, this set of 
algebraic equations for the modes a^, and therefore 
determining the apparent horizon as described by 

Eq. (nni). 



III. NUMERICAL RESULTS 

In this section, we present the numerical tests of 
our code as well the results about the dynamics of 
apparent horizons in RT spacetimes. We first need 
to specify the initial data K{u — 0,x) that deter- 
mine the initial values for the + 1 modes 6^(0) 
through 



2(lnA-(0,a:),P,) 



(21) 



We are going to consider two initial data in our nu- 
merical experiments. The first represents the ex- 
terior spacetirne of a homogeneous oblate spheroid 
described by 



K{0,x) = 



l + ^{a{Co) + lpiCo)P2{x)) 



(22) 



where and Bq are free parameters and a (Co) — 
arctan(l/Co), PiCo) = (1 + 3Co') arctan(l/Co) - Ko- 
There is a clear astrophysical motivation for such a 



family of initial data as pointed out in the works 
on the axisymmetric gravitational collapse of oblate 
gas spheroids satisfyiiig the Vlasov equation either 
in Newtonian theory [20|, as well in its relativistic 
generalization [2l|; and also connected with the ef- 
ficiency of the emission of gravitational waves [2^ . 
The second initial data family describe two initially 
boosted Schwarzschild black holes with opposite ve- 
locities V — tanh 770 [13i\ in which 



K{0,x) = ( 



\J cosh rjo — X sinh rjo 



Ao 



\J cosh 770 + X sinh 770 



(23) 



where A\ and Ai are arbitrary positive constants 
associated to the mass of each black hole. 

We first exhibit an important test for the spectral 
code used to integrate numerically the field equa- 
tions (l2|) and ([3]). In spite of non-stationary ana- 
lytical solutions of the field equations are not known 
(unless in the linear regime) , there exists a conserved 
quantity 



/o = 



K^{u, x)dx, 



(24) 



which is derived from the field equations and in- 
terpreted as the area of the fundamental 2-sphere 
spanned by {0,(j)). The conservation of Iq can be 
deduced after multiplying Eq. ^ by and inte- 
grating in the angular domain. Then, by specifying 
the initial data K{0, x), the initial value of Iq will be 
fixed and must be kept constant until the asymptotic 
state is achieved. In order to test if the numerically 
generated solution is accurate, we have evaluated the 
relative error between the numerical and exact val- 
ues of /o, (7 = |/o — /numer|/-foj whose rcsult is shown 
in Fig. 1 where we have included the influence of in- 
creasing truncation orders, or A = 7, 9, 13 (cf. Eq. 
(O). As it can be observed, the conservation of Iq is 
attained to about 10~^% accuracy for the smallest 
truncation order N = 7, and about 10~^^% accuracy 
for A^ — 13. Therefore, this result is a vivid proof 
of the accuracy of the numerical evolution scheme 
despite those other tests described in Ref. [Tsj . 

We now proceed with the numerical tests of the 
algorithm used to solve the apparent horizon equa- 
tion ([5]) at each time level u. Two steps will be 
needed. The first is to evaluate the time evolution 
of all modes bk (u) by integrating the dynamical sys- 
tem resulting from (|12p . In the second step these 
modes calculated at each u are inserted into the sys- 
tem of M + 1 algebraic equations derived from (20) 
and solve them to obtain the corresponding modes 
afc(u) that describe the apparent horizon through 




FIG. 1: Evolution of the relative error a between Jo 
and /numcr- We have considered the initial data that 
represent two initially boosted Schwarzschild black holes 
(Eq. ([231)) with Ai = 1.0, A2 = 0.2 and rjo = 0.3, as weU 
distinct truncation orders A'^ = 7, 9, 13 associated to the 
Galerkin expansion of K{u,x) given by Eq. ([9]). 



the function V{u,x) (cf. Eqs. HH) and In 
this way, the evolution of the apparent horizon is 
obtained until the stationary solution is attained. 
As a matter of fact, as an important piece of evi- 
dence of the accuracy of our numerical scheme, we 
have plotted in Fig. 2 the modulus of the residual 
equation (I16|) evaluated at the initial instant u = 
for both initial data families (|^^ and taking 
into account distinct values of the truncation orders 
N and M associated to the functions K{u, x) and 
V{u,x) (cf. Eqs. ([9]) and (fTS]) ). According to these 
plots the residual equation approach to zero as a con- 
sequence of increasing the truncation orders under 
consideration. 

A more enlightening experiment for depicting the 
convergence of the code is to exhibit the evolution of 
the L2 norm corresponding to the residual equation 
(dni) given by 



ResAH('u, x)^dx, 




(25) 



considering again distinct values of the truncation 
orders N and M. From Fig. 3 it can be seen that 
the norm evaluated at it = 0.4 decays exponentially 
if the truncation order N is increased, which demon- 
strates the expected geometric convergence typical 
of spectral methods. In Fig. 4 the full evolution 
of L2 is presented for increasing truncation orders 
M, N, and as expected it is noticed a rapid decreased 
of the norm until reaching to the value considered 
zero up to our numerical precision; also when the 
truncation order is increased, less time is necessary 




A'= 16 
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FIG. 2: Log-linear plots of the residual equation associ- 
ated to the apparent horizon equation (|16p showing the 
effect of increasing the truncation orders A'^ as indicated 
above and Af = A'' in each case. Notice the decrease 
of the residual equation as more terms are included in 
the Galerkin expansion. The initial data are the oblate 
spheroid given by (|22|) with ^0 = 0.1 and Bo = 1 (first 
graph), and the two initially boosted black holes (Eq. 
(1231)) in which Ai = 1.0, A2 = 0.1 and rjo = 0.3. 



to reach to that value. 

The evolution of the apparent horizon is illustrate 
by a sequence of polar plots of r = V{u, x) and de- 
picted in Fig. 5. We started at u = with the oblate 
spheroid initial data (j^^ and several plots are shown 
in subsequent instants until Uf = 500 where a circle 
is formed. Indeed, this is a consequence from the fact 
that the asymptotic state is the Schwarzschild con- 
figuration characterized hy K = constant according 
to Eq. ^ which implies also in V — constant. 

An interesting application of our code is to fol- 
low the behavior of the apparent horizon mass that 
is basically the amount of mass enclosed by the ap- 
parent horizon. It is worth of mentioning that the 
apparent horizon mass has thermodynamical prop- 
erties similar to those associated to black holes [1^| . 
In the case of RT spacetimes the past apparent hori- 
zon can only decrease in area, and therefore its mass 
decreases, contrary to the monotonic increase of the 
future apparent horizon area. The apparent horizon 
area Sah is evaluated through the following expres- 



N 

FIG. 3: The L2 norm of the residual equation (|25|l evalu- 
ated at u = 0.4 and for distinct truncation orders. Notice 
that the exponential decay of the norm with the increase 
of the truncation order is typical of spectral methods. 
Here we have used the two boosted black holes (|23|) with 
Ai = 1.0, A2 = 0.15 and -qo — 0.4 as the initial data. 
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FIG. 4: Behavior of the L2 norm of the apparent horizon 
equation (cf. Eq. (|25|l 'l for the oblate spheroid (|22|l with 
^0 = 0.1 and Bo = 0.3 taking into account distinct trun- 
cation orders A'^ and with Af = A'^ in each case. Again 
as far as A'' is increased more rapidly the norm tend to 
zero. 



Sah ~ '^^ J r^K^{u,x)dx = 

= %nml e-^'"^'''''^K^{u,x)dx, (26) 

where r = V{u, x) = 2m,Qe~'^^^"'^^ describes the ap- 
parent horizon (cf. Eq. (dH)), and the apparent 
horizon mass is expressed as 

In Fig. 6 we present the evolution of the apparent 
horizon mass and the Bondi mass [H, 



3 K 

2 

FIG. 5: Dynamics of all modes akiu) through a sequence 
of polar plots of 2mo exp(— _Fa(M, x)) starting from the 
oblate spheroid initial data {C,o = 0.1 and -Bo = 0.4) 
at ?io = as indicated, and for several times until the 
Schwarzschild final configuration (represented by a cir- 
cle) is formed at approximately u/ = 500. 



Mb = -Too / K^{u,x)dx, (28) 
2 J-i 

for the first family of initial data P^ . Accord- 
ing to Ref. [l2| the asymptotic configuration is 
the Schwarzschild black hole whose final mass as- 
sumes the value Mbh ~ rn^KQ, where Kq — 
lim„^oo K[u,x). This amount is smaller than the 
mass associated to the initial data since part of it 
is extracted by gravitational waves [2^ during the 
evolution of the spacetime, and consequently pro- 
ducing a monotonic decrease of the Bondi mass as 
shown in Fig. 6. Nonetheless, the decay of the ap- 
parent horizon mass is due to the decrease in area 
of the past apparent horizon as we have mentioned 
before. Notice also that both final values of the 
apparent horizon and Bondi masses coincide. As 
a matter of fact, this result is expected from the 
asymptotic solution of the apparent horizon equa- 
tion 14sympt = ^ttiqKq, and together with the ex- 
pression for the apparent horizon mass (|27p one can 
arrive at Mah — Mbh = moKQ. 

The final task is to consider the second family of 
initial data which represent the head-on colli- 
sion of two Schwarzschild black holes and a gener- 
alization of the initial data that describe the 
exterior of an inhomogeneous oblate spheroid [T^ . 
The common feature of both initial data is that 
the asymptotic configuration will be a boosted black 
hole [13, U^l described by 
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150 
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FIG. 6: Numerical evolution of the Bondi and appar- 
ent horizon masses in units of mo attesting their mono- 
tonic decay until the final configuration identified as the 
Schwarzschild black hole in which both values coincides. 
The initial data is the oblate spheroid with <^o = 0.1 and 
Bo = 0.4. 



by a if -transformation of the BMS group [2g| cor- 
responding to a boost along the z-axis, as given by 
Eq. ((29)) . This is the origin of the gap in Figs. 
7, which has the value Mb/Mah = cosh/^, where 
fi is the boost parameter of the if-transformation 
specified in Now the entropy of each final 

black hole, considered as a thermodynamical sys- 
tem in equilibrium, is defined as proportional to the 
area A of its event horizon and is invariant by a 
if-transformation as can be easily verified. This 
should be expected since an eventual possible def- 
inition of the BH entropy by a counting of its mi- 
croscopic states could not depend, in principle, on 
the state of motion of stationary black holes rel- 
ative to inertial frames at infinity. Therefore we 

have s(^A/4:nG^ = SMb/Tb = 6Mah/Tah so that 
this gap also defines the temperature transformation 
Tb Tb cosh/i between the two inertial rest frames 
of the BHs. 



IV. FINAL CONSIDERATIONS 



lim K{u,x) = , ^ ■ , , (29) 
u^oo cosh ^ + X smh /i 

where the values of Kq and the boost parameter /i 
are fixed by the numerical solution of the RT equa- 
tion (see Ref. ^1] for details), and the final Bondi 
mass is given by Eq. ([6]). It is worth mentioning that 
the imbalance in momentum of the initial gravita- 
tional wave distribution is responsible for the boost 
of the resulting black hole. In Figs. 7(a) and 7(b) 
we observe again the monotonic decay of both Mah 
and Mb with the retarded time u, but there is a gap 
between their asymptotic values. In order to under- 
stand the origin of this gap, we have noticed that 
according to the numerical experiments the asymp- 
totic solution of the apparent horizon equation ([8]) 

- 2 

is the same as the previous case, Vasympt — Smoi^To : 
in spite of i^asympt not being a constant (cf. Eq. 

Therefore, the final value of the apparent 
horizon mass can be evaluated from Eq. (|77|) (note 
that in this situation K = K{x)), and whose result 

— 3 

is Mah = moKo . In fact, this value is exactly the 
rest mass of the boosted black hole and consequently 
the gap observed in both graphs of Fig. 7 is due to 
the kinetic energy of the resulting black hole. 

The above results can be interpreted in the light 
of the so called First Law of Black Hole Thermo- 
dynamics [2^- [13]. The final mass configurations 
displayed in each of Figs. 7 can actually be inter- 
preted as two static black holes boosted with re- 
spect to the each other, namely, they are connected 



In this paper we have implemented and tested 
a numerical scheme based on a combination of 
Galerkin and pseudo-spectral methods to solve the 
apparent horizon equation in RT spacetimes. This a 
direct extension of the previous algorithm [l^ used 
to integrate the field equations ^ and The 
apparent horizon equation is reduced to a set of 
nonlinear algebraic equations for the modes ak and 
whose solution at each instant determines the ap- 
parent horizon described by Eq. ((S]). The applica- 
tions have consisted in solving the apparent horizon 
corresponding to initial data describing the exterior 
fields of oblate spheroids and the collision of two 
Schwarzschild black holes. 

We have performed numerical tests that strongly 
indicate the convergence and accuracy of the code. 
In our numerical experiments two initial data fam- 
ilies in RT spacetimes were considered: the first 
represents the gravitational field outside a oblate 
spheroid while the second two initially boosted 
Schwarzschild black holes with opposed velocities. 
We have confirmed that the Bondi mass Mb de- 
creases monotonically as the result of the mass 
extraction due to the gravitational waves, until a 
asymptotic value that coincides with the total mass 
of the resulting Schwarzschild black hole. The ap- 
parent horizon mass Mah also decreases with re- 
spect to u in face of the decreased in area which is 
expected in the case of past apparent horizon. By 
considering the first initial data, the asymptotic val- 
ues of Mb and Mah coincide to the value of the 
total mass of the resulting Schwarzschild black hole. 
In this case the apparent horizon mass is exactly the 
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mass enclosed by the event horizon. On the other 
hand, if we take into account the second initial data a 
gap between the asymptotic values of both masses is 
observed similarly as noticed by Chow and Lun . 
The origin of the gap is associated to the final config- 
uration identified as a boosted Schwarzschild black 
hole for which the Bondi mass is the total mass- 
energy content that includes the rest mass ttioKq 
plus the kinetic energy, whereas the final apparent 
horizon mass is the rest mass. 



Finally, in spite of RT spacetimes being the sim- 
plest asymptotically flat radiating geometries, they 
can be potentially used as simple but useful theo- 
retical laboratories to study relevant features of the 
bounded sources emitting gravitational waves (see 



Refs. [IJI, (131, UJ and 28]), and also to test 
new numerical schemes such the one we have im- 
plemented here. The natural step in our research is 
to examine the evolution of apparent horizons in the 
case of general RT spacetimes, and also in more re- 
alistic frameworks such as for spacetimes with Brill 
waves. 
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FIG. 7: Numerical decays of Bondi and apparent horizon 
masses in units of mo. Here the final configuration is a 
boosted black hole and both values of mass do not coin- 
cide asymptotically. The Bondi mass given by Eq. ((6) 
gather the rest and kinetic energy of the hole, whereas 
the apparent horizon mass Mah is a measure of the rest 
mass. We have considered in the first pannel a non- 
homogeneous oblate spheroid initial data given by 

K{0,x) = [l + 0.05 (^q(Co)-H|/3(Co)P2(x)) +0.05 exp(a;- 

„-|2 

0.3) , Co = 0.1, whereas in the second pannel the head- 
on collision of two boosted black holes with rjo = 0.05, 
Ai = 1.0 and A2 = 0.99. 
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